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ABSTRACT 

We consider the evolution of millisecond radio pulsars in binary systems 
with a main-sequence or evolved stellar companion. Evolution of non¬ 
accreting binary systems with “eclipsing” milisecond pulsars was described 
by Kluzniak, Czerny & Ray (1992) who predicted that systems like the one 
containing the Terzan 5 PSR 1744-24A will in the future become accreting 
low mass X-ray binaries (LMXBs), while PSR 1957-1-20 may evaporate its 
companion. The model presented in the current paper gives similar results 
for these two objects and allows to obtain diverse evolutionary tracks of 
millisecond pulsars with low mass companions (black widows). Our results 
suggest that the properties of many black widow systems can be explained 
by an ablation phase lasting a few hundred million years. Some of these 
sources may regain Roche lobe contact in a comparable time, and become 
LMXBs. 

Keywords: milisecond pulsar - redback - black widow - binary evolu¬ 
tion - ablation - LMXBs - gravitational waves 


1 INTRODUCTION 

Millisecond pulsars are probably intimately connected with LMXBs, as was real¬ 
ized soon after their discovery: it was suggested that millisecond pulsars have been 
spun up in LMXBs and will end their history in the radio pulsar phase (Radhakr- 
ishnan and Srinivasan, 1982; Alpar et ah, 1982). However, with the discovery of 
the eclipsing pulsars it was realized that some millisecond pulsars currently ablating 
their companions may re-enter the LMXB phase in a later epoch (Bisnovatyi-Kogan, 
1989; Ergma and Fedorova, 1991; Kluzniak et ah, 1992). Recent discoveries of many 
ablating binary systems have led to a rekindling of these ideas, and to the necessity 
of explaining the evolutionary status of these black widows and redbacks, as they 
are called (e.g., Roberts et ah, 2014). 
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We are presenting an evolutionary model describing a binary system composed of 
a pulsar and its stellar companion. The model includes effects like gravitational wave 
emission by the binary, ablation of the companion, and pulsar spindown. In general, 
part of the ablated matter may accrete onto the neutron star and another part may 
leave the system. The computed evolutionary tracks begin with the pulsar turn¬ 
on at the conclusion of the standard epoch of accretion in a semi-detached phase. 
Throughout most of the computed evolutionary history, the separation between the 
pulsar and the companion star is large enough for the latter to be below its Roche 
lobe. Therefore the only mechanism of mass loss considered in our model is ablation 
by the pulsar wind. 


2 MODEL DESCRIPTION 

The period of a binary system including a pulsar of mass M and its companion of 
mass m is 

27Tj^M + m) 

where J denotes total orbital angular momentum. The rate of change of the com¬ 
panion mass m is assumed to be proportional to the spin-down flux 
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where E is the energy loss of the pulsar primary owing to its spindown, d is the 
separation between the primary and the secondary, and a, b are model dependent 
exponents. In the simple model assumed in Kluzniak et al. (1992) a = b = 0. 
However, in Brookshaw and Tavani (1995) one may find a = 1/6 and b = —4/3. 
We will adopt the latter values. The change of mass of the primary is in principle 
connected with rh as M = —jdm. The coefficient (3 describes how much of the mass 
lost by the companion is accreted by the neutron star, and how much is lost from 
the binary in a wind, thus 0 < /3 < 1 with 0 corresponding to no accretion and 1 to 
no wind. We will take /3 = 0. 

The change of angular momentum [first term in Eq. (5)] is connected with two 
processes: emission of the gravitational waves (GW) and mass loss from system. 
We take the rate of angular momentum loss to gravitational waves to be described 
by (e.g., Shapiro and Teukolsky, 1983) 
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If we assume that specihc angular momentum carried away by a wind escaping from 
the system is j = aMJ/[m{m + M)], we have 
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Both Eqs. (3) and (4) conribute to the rate of change of the angular momentum: 
J = Jrh + Jqw- By differentiating Eq. (1) with respect to time we get the rate of 
change of the period 
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Equations (5), (6), (7), (8) constitute a system of first-order ordinary differential 
equations, which we proceed to solve with various assumptions and different initial 
conditions. 

In the simple case of no accretion onto the primary star, negligible companion 
mass, TO << M, and hence negligible gravitational wave emission, the equations 
reduce to (Kluzniak et ah, 1992) 
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and can be easily integrated, yielding 
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With suitable initial conditions the evolutionary paths on the P vs. to plot described 
by this equation can be made to pass through the current positions of some of the 
known pulsars, e.g. PSR 1957-20 (see the Appendix, Eig. 2). 

The source of the energy driving the ablation process is pulsar spindown. From 
the magnetic dipole formula (e.g., Shapiro and Teukolsky, 1983) we have 


E = - 


6c^ 


( 11 ) 


where B denotes the surface magnetic field near the pole, is the pulsar spin rate 
(the pulsar period being Pq = 2tt/^), Ris the pulsar radius and 9 denotes the angle 
between the magnetic and the rotation axes (for simplicity we take sin^ 9 = \ and 
i? = 10 km). On the other hand we have E = where / is the moment of inertia 
of the pulsar. These two equations provide 
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where Oq = 0(0) is the initial angular velocity of the pulsar and r = —O(0)/O(0) 
is the characteristic age of the pulsar (at time t=0). Eqs. (11), (12), are used in 
Eq. (6) to find to as a function of time. 
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When the secondary star is sufficiently close to the pulsar that it fills the Roche 
lobe, accretion through the inner Lagrangian point starts. This situation is not 
described by our model, although our tracks may bring the system to this point. 
When the radius of the companion is equal to Roche lobe radius, the relation 
between orbital period P and companion mass m is 

P = (13) 

where B « 0.462. The values of n and A correspond to the radius of the companion 
thorugh r = AmA. For degenerate stars like white dwarfs, n = —1/3, and for 
a hydrogen white dwarf A = 2.82 • 10‘^Mq^^ km (Shapiro and Teukolsky, 1983), 
while from Hamada and Salpeter (1961) one obtains A = 8.80 • IO^Mq ^ km for 
a helium white dwarf. One may also obtain this coefficient for a carbon white 
dwarf, which is A = 8.72 • 10 ^Mq' km (Hamada and Salpeter, 1961), it is almost 
indistinguishable from the helium one. Lines corresponding to Eq. (13) indicate 
where the evolutionary track may terminate in a Roche-lobe overflowing LMXB, 
depending on the companion type (Fig. 1). 


3 RESULTS OF NUMERICAL CALCULATIONS 

Using Mathematica, we solved numerically the system of four differential equations, 
Eqs. (5), ( 6 ), (7), ( 8 ) discussed above. We consider a model with no accretion 
(/3 = 0), we assume 7 = 2.5 • 10^ s^°/^ g“^/® (cf., Chen et ah, 2013) and, following 
Kluzniak et al. (1992), we take a = 0.86. For the initial point on the (m, P) plane 
we use one of two points on the track of Tauris and Savonije (1999), which describes 
the evolution of a LMXB with an evolved companion. For the PSR 1957-1-20 and 
B1744-74A (Terzan 5) tracks we use the starting point of Kluzniak et al. (1992), 
corresponding to the point at which magnetic braking is supposed to lose importance 
in the evolution of binaries with a main sequence companion. Current system 
parameters are taken from ATNF Database (2014); Manchester et al. (2005), and 
they can be found in Table 1, together with other data, for the six tracks which are 
presented in Fig. 1. 

Derived times of evolution are tgv — 7-10® y for PSR 1957-1-20 and tev — 5.5 • 10® y 
for Terzan 5. For PSR 1957-1-20, evolution is steady, whereas for Terzan 5 one can 
distinguish three stages of evolution. The first stage, when the evolution curve is 
nearly a straight line, lasts about 4.5 • 10® y. The second one, when the evolution 
path “turns downwards” on the P vs m plot, lasts 2.4 • 10® y. The last stage, 
when gravitational radiation is dominant, lasts 2.6 • 10® y. Objects with convex 
evolution curves evolve comparably fast: e.g. for J1807-2459A the evolution time is 
fev — 5.7 • 10® y. The values of tev in parentheses in Table 1 (for the Terzan 5 pulsar 
and J1023-I-0038) correspond to the time it will take for the system to regain the 
line of Roche-lobe contact starting from the present position. 
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Quantity 

PSR 1957+20 

B1744-74A/Terzan 5 

Initial 

Present 

Initial 

Present 

P[hr] 

2.9 

9.2 

3.0 

1.82 

m[M©] 

0.235 

0.022 

0.235 

0.087 

M[Mq] 

1.7 

1.7 

1.4 

1.4 

m[g/s] 

-2.0-103^ 

-7.6-1033 

-1.1 -1037 

-3.1 - 1033 

Po[ms] 

0.92 

1.60 

1.95 

11.56 

B[G] 

4.0- 10® 

1.3-109-6 

p[G- cm3] 

4.0 ■ 1026 

1.31027 

tev [y] 

6.72-106 

5.52-10® (5.90-10®) 

Quantity 

J1807-2459A 

J2241-5236 

Initial 

Present 

Initial 

Present 

P[hr] 

1.75 

1.71 

1.75 

3.50 

m[M0] 

0.115 

0.009 

0.115 

0.012 

M[Mq] 

1.4 

1.4 

1.4 

1.4 

m[g/s] 

-4.5 • 1036 

-2.6-1036 

-1.4-1036 

-1.1 - 1036 

Po[ms] 

2.91 

3.06 

1.98 

2.19 

P[G] 

2.9- 10® 

2.4-10® 

^[G-cm3] 

2.9 • 1026 

2.4 -1026 

tev [y] 

5.60-106 

7.88-10® 

Quantity 

J1311-3430 

J1023+0038 

Initial 

Present 

Present 

Predicted 

P[hr] 

1.75 

1.56 

4.73 

0.31 

m[M0] 

0.115 

0.008 

0.136 

0.061 

M[Mq] 

1.4 

1.4 

1.4 

1.4 

m[g/s] 

-3.6-1036 

-3.0-1036 

-2.2-1036 

LMXB 

Po[ms] 

2.48 

3.56 

1.67 

11.09 

P[G] 

2 .0- 10® 

7.9-10® 

p[G ■ cm^] 

2.0 -1026 

7.9 -1026 

tev [y] 

5.44-10® 

(1.8-1036) 


Table 1. System parameters 
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Figure 1. Large dots correspond to the present parameters of the observed pulsar systems 
for which the evolutionary tracks have been computed (dashed blue lines). Small dots are 
other objects taken from the ATNF Database (2014). Solid lines correspond to Roche lobe 
contact for a cold companion. The short-dashed red line, taken from Tauris and Savonije 
(1999), corresponds to LMXB evolution of a system with an evolved companion. Also 
shown (thin horizontal line segments) are the positions of two short-period LMXBs. The 
filled triangles mark plausible initial points of the evolutionary tracks. 


4 EVOLUTIONARY TRACKS 

Evolution of the system depends on the ratio between angular momentum losses 
caused by ablation and gravitational wave emission. There seem to be three types 
of tracks. 
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In the case where gravitational waves emission can be neglected (like in the PSR 
1957-20 system) the track is well described by the formula of Eq. (10). The system 
very nearly follows a straight line on a logP versus logm plot. The slope of this 
line depends only on the parameter a. The track may be deflected a little bit due 
to vestigial gravitational wave emission. 

Another possible track passes through the Terzan 5 pulsar B1744-24A. In the 
initial phase of system evolution the track is similar to the one described in the 
previous paragraph. The difference is that at a certain moment, owing to pulsar 
spindown, gravitational wave emission starts to dominate over ablation. If, from 
that point on, mass loss were neglected (i.e., the evolution were driven by GW 
emission alone), the track would be a vertical line on the logP-logm plot. In 
fact, a residual effect of ablation is still felt, and the track deviates slightly in the 
direction of lower companion mass (to the left in the figures). 

Neglecting mass loss from and mass transfer in the system (7 = 0 in Eq. [ 6 ]) one 
easily obtains the time elapsed in the evolution from binary period Pi to period P: 


T = 


5c^ 

20487r3GJi 


(p8/3^1/3 _ 


(14) 


where Ji is the initial angular momentum (corresponding to Pi). Time scales of 
evolution obtained from this equation are similar to the numerical values for the 
nearly vertical tracks in Fig. 1. 

Tracks similar to those described above were already obtained by Kluzniak et al. 
(1992). They cover situations in the limit where one of the effects, ablation or 
GW emission, dominates over the other along each major segment of the trajectory 
(although, as remarked above in Section 3, PSR B1744-74A spends most of its evo¬ 
lutionary time in transition between two such states). It seems that systems with 
an evolved very low mass companion (m < 0.04Mq) cannot evolve this way. For 
instance, obtaining a “Terzan-like” evolution track for these systems leads to evolu¬ 
tion time amounting to several dozens of billion years. A third type of evolutionary 
track seems to be required. 

We have found evolutionary tracks connecting the currently observed binary pa¬ 
rameters of the pulsars J2241-5236, J1807-2459A and JI3II-3430 with a plausible 
inital point and having reasonable time scales of evolution. These evolutionary 
tracks are characterized by angular momentum loss to both GW emission and ab¬ 
lation effects, and have a convex shape on a log P versus log m plot (Fig. I). Even¬ 
tually, the separation of the system components becomes large enough that GW 
emission loses importance, and the track becomes parallel to that of PSR 1957-1-20. 


5 DISCUSSION 

We have considered the evolution of millisecond radio pulsars with binary low-mass 
companions assuming simple formulae for the ablation rate of the companion by the 
pulsar wind. For the starting point of each evolutionary track that we considered we 
have taken a plausible moment of pulsar turn-on in an erstwhile LMXB, either along 
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the standard evolutionary curve familiar from discussion of cataclysmic variables 
and the period gap, i.e., a binary with a main-sequence companion (Paczyhski 
and Sienkiewicz, 1983), or along an evolutionary track with an evolved companion 
Tauris and Savonije (1999). Pulsar turn-on (or turn-off) in (potentially) accreting 
low mass binaries was discussed in Kluzniak et al. (1988). 

We have reproduced the results of Kluzniak et al. (1992) who performed a simi¬ 
lar study for the only two known eclipsing pulsars at the time (PSR 1957-1-20 and 
B1744-74A in Ter 5), and found that there are periods of their evolutionary history 
in the ablation phase when either one or the other of two major angular momen¬ 
tum loss mechanisms dominates (mass loss from the system or GW emission). We 
note that evolutionary tracks that we now find based on the Brookshaw and Tavani 
(1995) evaporation formula, Eq. (2), imply shorter initial pulsar periods than pre¬ 
viously obtained, this can be seen from a comparison of the entries in Table 1 with 
the description of tracks (b) and (e) in the Appendix, Fig. 2. 

We find that we are able to reproduce the current positions of typical millisecond 
radio pulsars with a low mass binary companion, typically this involves an ablation 
phase lasting several hundred years. However, we find that for the majority of the 
black widow pulsars known today the relative importance of the two considered 
angular momentum loss mechanisms is comparable in their evolutionary history, 
i.e., unlike in the case of PSR 1957-1-20 and B1744-74A, neither GW emission nor 
mass loss domiantes the other over major portions of the evolutionary track in the 
period-mass diagram (Fig. 1). 

We conhrm the conclusion of Kluzniak et al. (1992), who predicted that some ms 
pulsars may become accreting LMXBs at the end of their evolution. Two of the 
tracks presented in this paper end very close to the line of Roche-lobe contact, in 
the current position of PSR J1807-2459A and PSR J1311-3430. These two pulsars 
seem to be close to the end of a 5 x 10® y ablation phase. 

We note that detailed binary evolutionary calculations, which included an abla¬ 
tion phase similar to the model considered here were presented recently in Ghen 
et al. (2013). 
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6 APPENDIX 

For ease of reference, we reproduce Figure 1 and its caption from the pre-arXiv 
contribution of Kluzniak, Czerny and Ray (1992). 



Figure 2. ‘Figure 1. Possible evolutionary tracks of systems with “evaporative” mass loss. 
The decimal logarithm of the orbital period in hours is plotted versus the decimal logarithm 
of the mass of the companion in units of Solar mass. Likely location of the eclipsing pulsars 
(filed circles) as well as possible positions of the X-ray binaries 4U 1916-05 and 4U 1820-30 
are also indicated (dash-dot-dash lines). The thick straight line segments correspond to 
systems with a main-sequence or a cold degenerate dwarf companion in Roche-lobe contact. 
According to the standard theory of their evolution, cataclysmic variables follow the thin 
curve (in the direction of decreasing companion mass, m). When this theory is applied to 
canonical LMXBs, the dotted tracks ensue, see Section 5 for details. The lines (a) through 
(e) differ only in the properties of the pulsar ablating its companion: in the strength of the 
magnetic dipole moment and in the initial value, Pq, of the rotational period of the neutron 
star. The values of Pq and log(i3/Gauss), where B = fix 10“^® cm“®, are respectively (a) 
5.0 ms, 9.5; (b) 3.4 ms, 8.9; (c) 2.0 ms, 9.0; (d) 2.0 ms, 8.6; (e) 1.25 ms, 8.1. We assumed 
that 10% of the energy flux impinging on the companion is converted into kinetic energy 
of the evaporative plume, and we took /3 = 0.86.’ 


N.B. The parameter “/3” in the quoted caption corresponds to our a. 





